Euler-Lagrange correspondence of generalized Burgers cellular 

automaton 



Junta Matsukidaira and Katsuhiro Nishinari 
Department of Applied Mathematics and Informatics, 
Ryukoku University, 
Shiga 520-2194, Japan 

(Dated: February 8, 2008) 

Abstract 

Recently, we have proposed a Euler-Lagrange transformation for cellular automata(CA) by de- 
veloping new transformation formulas. Applying this method to the Burgers CA(BCA), we have 
succeeded in obtaining the Lagrange representation of the BCA. In this paper, we apply this 
method to multi-value generalized Burgers CA(GBCA) which include the Fukui-Ishibashi model 
and the quick-start model associated with traffic flow. As a result, we have succeeded in clarifying 
the Euler-Lagrange correspondence of these models. It turns out, moreover that the GBCA can 
naturally be considered as a simple model of a multi-lane traffic flow. 

PACS numbers: 05.45.Yv, 05.50.+q, 89.75.-k, 89.40. +k 
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INTRODUCTION 



Cellular automata(CA) are simple objects that can describe not only physical systems 
but also complex systems such as chemical, biological and social systems. Since Wolfram has 
introduced the elementary CA (ECA) in one-dimension as a model of complex systemjlj], 



lots of studies have been done in these fie 
Anions: them, conservative CA 



ongtJ 

SB, 



ds. 



3, S, E|, including soliton systems P, [f| and traffic-flow 
models[7|, |8(, have attracted much attention because of their analytical properties which can 
be treated by various mathematical methods. 

In recent years, the authors have proposed and studied the so-called ultra-discrete method 
by which we obtain corresponding CA from difference equations^, 0,0,0, HQ. Various 
soliton equations have been transformed into CA by this method and the solutions of the CA 
have also been obtained exactly without losing the mathematical properties of corresponding 
difference equations. In this method, max-function plays an important role through the 
ultra-discrete formula, 

lim^elog ^exp ^— j + exp ^ — jj = max(v4, B). (1) 

On the other hand, it is known that several conservative CA allow two different represen- 
tation, Euler representation and Lagrange representation [15]. In the Euler representation, 
particles are observed at a certain fixed point in space as a dependent (field) variable, while 
in the Lagrange representation, we trace each particle and follow the trajectory of it. Thus 
a dependent variable represents the position of each particle in the Lagrange representation. 
In the previous paper |l6||, we have proposed a Euler-Lagrange transformation for CA by de- 
veloping new explicit transformation formulas. We started from Burgers CA(BCA), which 
is the Euler representation of rule-184 CA, 

Vf X = U) + min(^_ 1 , 1 - U*) - min(C/*, 1 - (2) 

where Uj denote the number of particles at the site j and time t and introducing the variable 
Shy 

k=— oo 
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or Uj = — Sj_ 1: and putting 

N-l 

fij = E^'-^ ( 4 ) 

i=0 

where H(x) is the step function defined by H(x) — 1 if x > and = otherwise, we 

obtain Lagrange representation of rule-184, 

x\ +1 = min(x* + l,x* +1 -l) 

= x* + min(l,x* +1 -x\- 1), (5) 

where is the Lagrange variable that represents the position of the i-th. particle at time t. 

In deriving the Lagrange representation from Euler representation, following two formulae 
for max function and step function play important roles. 

One is 

n / n n \ 

^2 H {j - min(a fc , b k )) = max I ^ H U ~ a k), E H ^ ~ b ^ ) ' ( 6 ) 

k=l \k=l k=l J 

where we assume that a± < a<i < ■ ■ ■ < a n and b\ < hi < ■ ■ ■ < b n . This formula expresses 
the commutability of max function and step function. Another formula is 



max H(j - a*) - m, = £ H(j - a< + J, (7) 

where we assume that a* = oo if j is larger than the number of particles. 

In the Euler-Lagrange transformation, max-function again plays an important role as 
in the ultra-discrete method. In addition, the step function significantly contributes the 
transformation as well. 



GENERALIZED BURGERS CELLULAR AUTOMATON 

In this paper, we present a new model which is a multi-value (multi-lane) and multi- 
neighbor (high-speed, long-perspective) extension of BCA. The equation in the Euler repre- 
sentation is expressed by 

(v-i p-i \ /v-i P-l \ 

E u U-k' E( L - u hk) - min E U U> E( L - ^+1+*) , (8) 
k=0 k=0 J \k=0 k=0 J 
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where Uj G {0, • • • , L} denotes the number of particles at the site j and time t. The 
parameter L represents the maximum capacity of a cell, V represents the maximum speed 
of particles and P represents perspective of particles, i.e., the maximum number of particles 
that a particle can see in its front. We call this model generalized Burgers cellular automaton 
(GBCA) in this paper. It will become clear later that GBCA include Euler representation 



and quick-start model 



rJQ, 



which are both known as traffic 



of Fukui-Ishibashi model 
models. 

Let us derive the Lagrange representation of GBCA according to the method which we 
have proposed in the previous paper [3]. First we introduce the variable S by 

S}= E Ut (9) 

k=— oo 

or Uj — — S t j_ 1 where S*j is the total number of the particles from — oo to j-th site, which 
is assumed to be finite. Rewriting © in S, we obtain 

Sj +1 =rna X (Sj_ v ,Sj +P -LP). (10) 

Here, we put 

N-l 



i=0 



where H(x) is the step function and N is the total number of particles on the cells. x\ is the 
Lagrange variable that represents the position of the i-th particle at time t and the relation 
Xq < x\ < • • • < xJv.x holds. 

Using (fTT|) to replace S in (fTU|) by H, we have 

N-l /N-l n-i \ 

E h v - = max E h v y )' E h v -A + p)-Lp\. (12) 

i=0 V i=0 i=0 / 

By using (JHJ) and (JTj), we have 

E H(j - x* +1 ) = max ( E H{j - x\ - V), E H(j - x\ +LP + P) 

i \ i i / 

= E^O- min & + y ' x ^p ~ P)) , (13) 

i 

Comparing both sides, we finally obtain 

x* +1 = min(x* + V, x\ +LP - P) 

= x\ + min(V, x\ +LP -x\-P) (14) 
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This is the Lagrange representation of GBCA. If we put L = 1,P = 1, ()14|) becomes the 
Lagrange representation of the Fukui-Ishibashi model|j|, 



x 



t+i 



x\ + m.m{V,x"-x\-l). (15) 



Therefore we know the Euler representation of Fukui-Ishibashi model is 

E/J+ 1 - V) = min (j2 U U-k> 1 - Uj ) - min ( E U U> 1 - U UA ■ ( 16 ) 
\fc=o / \fe=o / 

If we put L = 1, V = 1, we obtain the Lagrange representation of the quick-start mo del 1 171. 

m 

z* +1 = 4 + min(l,x* +P -x*-P). (17) 
Corresponding Euler representation of the quick-start model is 

U] +1 -V* = min \IJ)_^ - U* +k )\ - min (u], ^(1 - ) . (18) 

V fc=0 / V fc=0 / 

The expression in the case of P = 2 has already been obtained in 18]. Thus (HHD is a 
generalization of the previous result to the arbitrary P. 



MULTI- VALUE (MULTI-LANE) CA 

In the previous section, we show that (JBJ in the case of L = 1 includes the Euler rep- 
resentation of the Fukui-Ishibashi model and the quick-start model, and clarify the Euler- 
Lagrange correspondence of these models. 

In this section, we show that (jSJ) in the case of L > 1 can be interpreted as an multi- value 
(multi-lane) extension of the Fukui-Ishibashi model or the quick-start model. 

Let us consider (jHJ) in the case of P = 1 for simplicity. The Euler representation becomes 

f/f 1 -U) = min (fl U U-k> L ~ U j)- min ( E U U> L ~ Vi+i) . ( 19 ) 
\k=o J \k=o J 

and the corresponding Lagrange representation becomes 



x 



* +1 = x* + minO/,x* +L -x*-l). (20) 



We show in the following that (|19J1 is a (L + l)-valued CA and (|2(J|) describes the movement 
of particle in L-lane, which means ()19|) and ()20j) are considered as a (L + 1)- valued (L-lane) 
extension of Fukui-Ishibashi model. 
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Assume that L > and < C/| < L for any j at a certain £. Then, relations 



'V-l 



mm 



mm 



mm 



mm 



£?J-i-*> L-U] 



K k=0 
'V-l 



E^U 



> o 

> 



,fc=0 

<v-i 



v-i 



k=0 / 
V-l \ 



mm 



>fe=o 



E vj-i-k, L )< L 
fc=-i / 

,fc=i 

v-i N 
< min ( ^ L — Uj 



(21) 



mm 



holds. Therefore, < E/| < i holds for any j. This means the equation (|19|) under the 
above condition is equivalent to a CA with a value set {0, 1, . . . , L}. 

In the case of L = 2 and V — 2, the time evolution of (JT§j) is illustrated by the FIG. 
1, where the number of the particles in the jth cell is equals to the value of Uj. We can 
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FIG. 1: time evolution of the equation (|19|) 
consider this CclSG clS db 3-valued CA in the Euler representation. 
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Next, we introduce an ordering line on cells as shown in FIG. 2 . We number particles 
along this line as shown in FIG. 3, where i denotes suffix of variable x\ in (|20p. For example, 
particles of the initial condition in FIG. 1 are numbered as FIG. 4. 



FIG. 2: ordering line 
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FIG. 3: numbered particles along ordering line 
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6) : S 



FIG. 4: numbered particles at t = of FIG.l 

If we take FIG. 4 as a initial condition for (}2T)j) and move particles according to (f2*U|) . 
we obtain FIG. 5. The configurations of particles in FIG. 5 is identical with FIG. 1 if 
we neglect numbers assigned to particles. Note here that although (|2(J|) is a form of a car 
following model, the exclusion principle in one cell, which means Xi ^ Xj holds for any i,j 
, does not hold any more. Actually we see x\ = x\ surely holds in FIG 5. Without the 
exclusion principle, ()20|) works successfully together with (j!9|) . and describe the movement 
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of the particles. Note that there are no overtaking particles in the course of time on the 
ordering line. As the the number of particle in one cell remains less than 2, the model is 
considered as a car following model in 2-lane in this case. 
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FIG. 5: time evolution of the equation (|2U)) 

Generally, as (JT§j) is (L+l)-valued CA and (|20|) is a car-following form of equation like the 
Fukui-Ishibashi model, we can consider (fT^j) and (j2*U|l as (L + l)-valued (L-lane) extension 
of Fukui-Ishibashi model. 

In the case of the quick-start model, we can similarly construct multi-value (multi-lane) 
extension of the model although we don't discuss in detail in this paper. 

FUNDAMENTAL DIAGRAMS 

In this section, we discuss about fundamental diagrams of GBCA. In the following, we 
consider a periodic boundary condition. GBCA can be expressed in a conservation form 
such as 

A t E/J + A 3 q] = 0, (22) 

where At and Aj are the forward difference operator with respect to the indicated variable 
and qj represents a flow. Average density p and average flow Q t over all sites are defined by 
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p 



(23) 




where K is the number of sites in a period. Since GBCA is in a conservation form, average 
density does not depend on time and we can use p without a superscript t. 

As we have shown, GBCA can be expressed by (JTUJ). As (jlUj) can be transformed into 
ultra-discrete diffusion type equation, asymptotic behavior of GBCA is similar to that of 
BCA. Since a pattern of the maximum flow is given by • • • 11 ■ • ■ 11 PO ■••00/ • • , the density 



p 



v 



and the flow at the phase transition from free to congested state become 



P 



and 



PV 



P + V 



P + V 



respectively la ]. Thus the fundamental diagram is given by FIG. 6. 
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FIG. 6: 
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CONCLUDING DISCUSSIONS 



In this paper, we have proposed GBCA which is a multi-value (multi-lane) and multi- 
neighbor (high-speed, long-perspective) extension of BCA. We have derived the Lagrange 
representation of GBCA, which includes the Fukui-Ishibashi model and the quick-start model 
when L = 1. By using the Euler-Lagrange correspondence for GBCA, we have also obtained 
the Euler representation of Fukui-Ishibashi model. Furthermore, we have obtained multi- 
value (multi-lane) extension of the Fukui-Ishibashi model and the quick-start model, and 
shown the multi-lane interpretation by introducing ordering line on cells. It is worth men- 
tioning again here that the success of the Euler-Lagrange transformation in the multi-lane 
case is due to the relaxation of the exclusion principle in (|2Uj) . 

We have applied our method to more general cases than the case in the previous papers 
and it works successfully. There are more several conservative CA model which allows 
both Euler and Lagrange representation. For examples, soliton systems called box and ball 
system j^] are known to allow two different representation. Now we are progressing on finding 
the Euler-Lagrange correspondence for soliton systems and shall report them in the future. 

We believe that establishing the general Euler-Lagrange correspondence of CA will make 
a new development of the ultra-discrete method and studies of CA. 
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